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Generalized density dependence in Skyrme effective interactions is investigated to get 
>- ; forces valid beyond the mean field approximation. Preliminary results are presented for 
infinite symmetric and asymmetric nuclear matter up to pure neutron matter. 

O 

g i 1. INTRODUCTION 

q ■ It is commonly accepted that it does exist a relation between compressibility, effective 
mass and the density dependence of a given effective force. Studies with Skyrme forces [ 
P have shown that the incompressibility and the effective mass cannot be chosen 
independently once the analytical form of the (single) density-dependent term has been 
chosen. This has led to the p 1//6 density dependence in Skyrme forces like SkM* which 
allows a value of around 220 MeV close to that extracted from the experimental 
breathing mode analyses [0 El and an effective mass m* around 0.7m simultaneously. 
The present status about compressibility has been summarized by G. Colo during this 
conference [Hj. 

Recently, the density dependence of phenomenological effective interactions, such as 
Skyrme or Gogny forces, has been revisited in the context of beyond mean field cal- 
culations [ EI- Indeed, while a dependence of the interaction on the density is well 
established for calculations at the mean field level [E], no strongly motivated prescription 
exists when several mean fields are mixed as in the Generator Coordinate Method (GCM) 
and the Projected Mean Field Method (PMFM). First, an extension of the Goldstone- 
Brueckner theory has motivated the GCM and the PMFM from a perturbative point of 
view for the first time [[£]■ In this extended context, a generalized Brueckner G matrix 
summing particle-particle ladders has been defined and may be used as a reference from 
which phenomenological interactions in GCM or PMFM calculations should be approxi- 
mated. It is possible to simplify this in-medium interaction to extend the validity of the 
Skyrme force [|B1IH| and to identify approximately the density dependence originating from 
Brueckner correlations in the context of mixed nonorthogonal vacua [Ej. The renormal- 
ization of three-body forces has also been advocated as an important source of density 
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dependence in phenomenological interactions. Indeed, their renormalization through a 
density-dependent two-body interaction has been proved to be justified at the mean field 
level in some particular cases [ E] • As for the resummation of Brueckner correlations, 
the corresponding density turns out to be the one of the mean field which is calculated. 
Consequently, a single density-dependent term has often been used in phenomenological 
forces, e.g. Gogny [ 10J and Skyrme [ [S| effective interactions. On the other hand, the 
density dependence renormalizing three-body force effects has been shown to be different 
from the one taking care of Brueckner correlations when going beyond the mean field 
approximation [6 . In fact, the use of two different, and theoretically motivated, density- 
dependent terms happens to have non negligible effects on collective spectra [111). Based 
on this analysis, it is legitimate to redefine the Skyrme interaction at the mean field level. 
Including two density- dependent terms, each related to a given physical origin, will allow 
a proper extrapolation of the interaction to GCM and PMFM calculations. 

A first attempt is presented in this paper, where we concentrate on the practical advan- 
tages of having two different density-dependent terms. Relating them to the resummation 
of Brueckner correlations and to the renormalization of three-body forces is the aim of 
our on-going work. In this way, this paper is concerned by a systematic investigation 
of infinite matter, symmetric and asymmetric up to pure neutron matter, using multi 
density-dependent terms in Skyrme forces. Our analysis is focused on p 1//3 terms which 
are connected with a kp expansion of the Brueckner G matrix. 

2. STANDARD SKYRME EFFECTIVE INTERACTIONS 

When using Skyrme effective forces with a standard density dependence p a , where p is 
the total density, the adjustable parameter a is strongly correlated to the incompressibility 
and the effective mass in nuclear matter. This is illustrated on Figure Q for different 
standard Skyrme forces with various powers a. In all cases, the different parameters of 
the force have been adjusted in order to obtain the saturation density po = 0.16 fm~ 3 
and the energy per particle E/A = —16 MeV in infinite symmetric nuclear matter. Once 
these two properties are fixed, the relation between m*/m and is entirely determined 
by a. This well-known feature, discussed in [|T], yields to the conclusion that a value of 
a around 1, like in SIII, does not allow to reach a correct compression modulus [U\ I12j. 
Only values of a ranging from 1/6 to 1/3 allow for an acceptable set {m*/m, i^oo}- 

It is enlightening to study a bit further the relation between m*/m and in the case 
of the standard parameterizations of Skyrme forces. Introducing as in [[J the notation 
O s = 3ti + (5 + 4x 2 )t2, one has 

m* ( 1 m _ \ _1 
and 

3tf /3tt 2 \ 2/3 2/3 3 n /3vr 2 \ 2/3 5 . 9 a+l 
K ~ = -^{-) p ° + 8 Qs {— ) P ° + 16 a(a + 1) ^ ' (2) 
The energy per particle in infinite nuclear matter at saturation density is given by: 

E 3h 2 /3tt 2 \ 2/3 2/3 3 3 ^ /3tt 2 \ 2/3 5/3 1 a+1 
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Figure 1. Correlation between incompressibility, effective mass and density dependence 
of standard parameterizations of Skyrme effective interactions. 



The saturation density being determined by the condition on the pressure -P(po) = 
which gives (if po 7^ 0) 

* (^f> + f** + le. (^)V + !«.«, + 1K « = o. (4, 

The three equations (J2J), © and (@J) form a system with three unknown quantities to, t 3 
and S . Solving this system for a given set (p , E/A, K^) provides the coefficients to, t 3 
and S , this latter determining uniquely the effective mass m*/m through Eq. ((TJ. 
Actually, the previous conclusion is incorrect in one special case. To see that, let us set 



F(e s ,t 3 ) = - 



e. 



3tt 2 \ 2/3 5 
2 J 3 s 



(5) 



It is trivial to check that if and only if a = 2/3 we have 



3ft 2 /3vr 2 \ 2/3 2 /3 3 3 , 5/3 E 

S 2 /3tt 2 \ 2/3 2/3 3 , 5/3 

In this particular case, there are three equations and two unknown quantities, t and 
^(©s;^)) an d, in general, this system has no solution. However, we can use the last two 
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equations to determine to and F(0 s ,i 3 ) and then calculate the corresponding value of 
Koq. Because of the definition (J3J), S can be freely chosen. This situation is illustrated 
on Figure n by the vertical dashed line, when a = 2/3, the value of is fixed and the 
effective mass can be freely chosen. However, this case is not interesting since it does not 
correspond to a realistic value of K^. 

3. MODIFIED DENSITY DEPENDENCE IN SKYRME FORCES 

Following several previous attempts [ H31 Clj and adding a density dependence on 
each term of the force, one can investigate the properties of generalized Skyrme effective 
interactions written as: 

V(r h r 2 ) = ^oiP i/3 (l+zovP«x)<Kr) 

i 

+ ^t 2 ,p^ 3 (l + a; 2j P CT )P'-5(r)P (7) 

j 

P' x <5(r)P 

k 



5(r)P' 2 + P 2 5(r) 



+ W 0kP k/3 a ■ 



where p = p(R) is the total density and see ref. [Q] for the other notations. Standard 
parameterizations discussed in Section ^correspond to j = k = and % = 0, 3 in SIII and 
i = 0, 1/2 in SkM* and SLy4 for instance. 

This generalized parameterization increases the number of parameters in the force, 
making it more phenomenological and reducing its predictive power. Besides, the fit- 
ting procedure may become intractable if most of the parameters are free. In this first 
exploratory work, we limit our study to the case with t±j = t 2 j = except for j = 
(hereafter we omit the second index) and three non zero values for £ i> namely i\, % 2 and 
i 3 out of the set {0,1,2,3}. The corresponding forces will be labeled by «2, ^3]- The 
spin-orbit term will not be discussed here since it has no contribution in infinite nuclear 
matter. 

The problems faced with the Skyrme forces with one density-dependent term are cured 
with this family of parameterizations. Following the same procedure to adjust the coeffi- 
cients, we can fix the saturation properties of infinite symmetric nuclear matter p , E/A, 
and m*/m independently and determine the coefficients of the force. 

However this procedure only guarantees to have realistic properties of nuclear matter 
in the vicinity of the saturation point. If we want a force which gives an equation of 
state in agreement with ab initio calculations at low density as well as high density 
(up to ~ 3p ) one should adopt a different fitting procedure. First, one can fit a set of 
points [pi, E/A(pi)] i=1 N which samples a realistic equation of state of nuclear matter. 
Then, the parameters of the force will provide the characteristic properties of nuclear 
matter (p , E/A, and m*/m) using equations (JTJ) to (jlj). In this kind of approach, 
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the presence of a density-dependent term with the power 2/3 in the force is particularly 
interesting. Indeed, the coefficient Q s is only constrained by the function F(Q S , t 02 ) whose 
value comes out of the fit. So, by choosing correctly the parameter t 02 one can obtain any 
desired value of the effective mass. 

4. RESULTS USING THE GENERALIZED SKYRME FORCES 

In this section, we present some results using the generalized forces [0, 1, 2], [0, 1, 3] and 
[0, 2, 3] for symmetric nuclear matter as well as pure neutron matter. For this purpose, 
we need to determine the coefficients toi, xoi and the two combinations of coefficients 
6 S and Q v = ti(2 + x±) + t 2 {2 + x 2 ), this latter expression is related to the isovector 
enhanced factor k [ occurring in the Thomas- Reiche-Kuhn sum rule of the isovector 
giant dipole resonance. The quantities E/A = —16 MeV, po = 0.16 fm~ 3 , Koo = 230 MeV, 
m*/m = 0.8 and k = 0.5 allow the determination of the parameters t^, Q s and 0„. The 
coefficient x i are determined by fitting the neutron matter equation of state provided by 
Akmal et al. [E3- The coefficients obtained accordingly are summarized in Tabled 



Table 1 

Coefficients t 0i and Xoi of the forces used in this section. See the text for the remaining 
parameters. 



Force 


^00 ^00 


toi Xoi 




^03 ^03 


[0,1,2] 
[0,1,3] 
[0,2,3] 


-1855.38 0.4006 
-1807.41 0.3291 
-1431.36 0.2828 


2343.12 0.7589 
2078.04 0.4094 


-488.28 3.6176 
3827.79 0.3814 


-299.81 3.4464 
-2650.09 0.7640 



Figure El shows the two EOS (energy per particle as a function of equilibrium density) 
for symmetric nuclear matter and pure neutron matter compared with the EOS of Akmal 
et al. [US]. The results are quite reasonable and rather similar for densities from to ~ 3po- 
The most salient feature of the parameterizations is that the coefficient of the term with 
the highest density dependence is always negative, leading to the pathological property 
for high density E/A — ► — oo. Nevertheless, as can be expected from the figure or check 
numerically, this problem shows up only at extremely large and unphysical densities and, 
thus, is not relevant for the applications to atomic nuclei and neutron stars. 

Another important property which can be extracted from this study is the symmetry 
energy aj. Its expression, obtained by a straightforward calculation [H], depends upon the 
coefficients t^ and xo«. On Figure El the symmetry energy is plotted as a function of the 
density and compared with the variational results of Lagaris et al. [El and the schematic 
parameterizations of Bao- An Li [HI)]. The three forces [0,1,2], [0,1,3] and [0,2,3] give results 
in good agreement with the variational approach, especially for p < 0.3 fm~ 3 , even though 
this quantity has not been directly fitted. For p = po, the three forces give exactly the 
same symmetry energy, a/ ~ 32.5 MeV, which is close to the commonly accepted value. 
This result is very encouraging and reinforce the idea that the two density-dependent 
terms in the Skyrme force open pertinent degrees of freedom. 
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Figure 2. Energy per particle in infi- 
nite matter function of equilibrium 
density. SNM: symmetric nuclear matter; 
PNM: pure neutron matter. Labels [0,1,2], 
[0,1,3], [0,2,3] refers to the density depen- 
dence of the forces: see text. Open dia- 
monds and triangles: EOS of Akmal et al. [ 

Eg. 



Figure 3. Symmetry energy in infinite 
matter as a function of the density. Li(a) 
and (b): phenomenological parameteriza- 
tions by Bao-An Li [fTtj|. see text. Open di- 
amonds: variational calculation of Lagaris 
et al. [Ej. Inset: symmetry energy in the 
vicinity of equilibrium density of symmet- 
ric nuclear matter. 



5. CONCLUSION 

We have shown that the choice of a standard Skyrme effective force with a modified den- 
sity dependence based on two terms enables to choose independently the incompressibility 
and the isoscalar effective mass. The gross properties of nuclear matter investigated here 
are quite reasonable: incompressibility (~ 230 MeV), isoscalar effective mass (~ 0.8m) 
and symmetry energy (~ 32.5 MeV). With this generalized density dependence we have 
been able to construct new Skyrme like parameterizations without collapse at relevant 
densities and which exhibit reasonable equation of state (EOS) for symmetric nuclear 
matter as well as for pure neutron matter compared to the recent realistic variational 
EOS of Akmal et al. [US]. 

Other choices can be explored for the density dependence provided they are physically 
well grounded. The choice p 2//3 for one of the density-dependent term lets total freedom for 
the effective mass. This feature is extremely interesting from the perspective of developing 
accurate and predictive force since it gives a control on the density of state around the 
Fermi energy where the correlations beyond the Hartree-Fock approximation can develop. 
Furthermore it implies only two additional parameters so that their total number remains 
quite small. 
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